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Status of the computation of /g^, £ and g 

Damir Bee ire vie 

Laboratoire de Physique Theorique, Universite Paris Sud, Centre d'Orsay, F-91405 Orsay-Cedex, France 

Current status of the computation of the neutral S-meson mixing amplitudes, with particular attention to the heavy-light meson decay 
constants, is reviewed. The values for these quantities, as well as for the coupling of the pion to the lowest doublet of heavy-light 
mesons, are given. 



Decay constants of heavy-light mesons enter in the most 
direct way the standard analyses of the CKM unitarity tri- 
angle [[l]|2|, namely through the frequencies of mixing in 
the neutral Z?-meson systems: 

Am d = C B m Bil flB Bd A 2 A 6 [(1 - p) 2 + ff] , 



Am d m Bll flfiih 2 \ n -,2 ^ -21 



CD 
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which, when combined with experimental values for the 
mass differences, constrain the vertex of the CKM trian- 
gle, through the circle in the p — fj plane. Pseudoscalar 
(and vector) meson decay constants, , /£> , are also 
important in the studies of the corrections to the factoriza- 
tion approximation in non-leptonic Z?-decay modes. These 
constants are also important ingredients in the research of 
low energy physics effects coming from physics beyond 
Standard Model [|5J. In what follows, I will focus on 
the computation of the pseudoscalar decay constants, al- 
though most of the discussion is equally applicable to the 
vector ones. I will briefly go through the novelties related 
to the bag-parameters B Bd/g , and close this short review by 
discussing the value of the phenomenologically important 
coupling of the pion to the lowest-lying doublet of heavy- 
light mesons. 

1 Strategies and difficulties to compute 
the decay constants 

The decay constant fu of the heavy-light pseudoscalar me- 
son H is defined through 



{0\QYvY5q\H(p)) = ip ll fH, 



(2) 



where Q is either c- or b-quark, and q is s- or u/d-qnark. 1 
Although the matrix element is the simplest one, the 

1 The mass difference between mj — m u quarks cannot be resolved by any 
of the methods to compute the non-perturbative effects. 



fact that one quark is heavy and the other is light makes its 
computation extremely difficult (in spite of the simplifica- 
tions stemming from the heavy quark and chiral symme- 
tries). 

In various quark models the decay constant is either a 
parameter of the model, or its value depends crucially 
on the specific choice of the model's parameters. This 
situation is clearly unacceptable as we seek a "model- 
independent" estimate (an ab initio determination), solely 
relying on the underlying theory, on QCD. A step in this 
direction was made by employing the duality sum rules, 
both in QCD and in heavy quark effective theory (HQET). 2 
The non-perturbative contributions in this approach are 
parametrized by the power corrections, the coefficients of 
which are the vacuum condensates. Apart from the chiral 
condensate, the essential non-perturbative input in evalu- 
ating fs comes from the gluon and the mixed quark-gluon 
condensates. Those two are either not well defined or their 
values are poorly known. For example, the gluon conden- 
sate, as estimated from the comparison of the sum rules 
with the corresponding experimental information on the 
p (n) resonances [ 5 1 and charmonium resonances [ 6|, differ 
by a factor of roughly 8: 



°^G* V ) = (0.074: 



0.023) GeV 



, (0.009 ± 0.007) GeV 4 



Besides, it is often unclear where the onset of the quark- 
hadron duality in the sum rule takes place (parametrized 
by the threshold parameter so). In short, the benefit of the 
method is that it made a step towards the first principle 
QCD calculations of the hadronic quantities; the drawback 
is the existence of too many parameters whose values are 
loosely constrained, thus prohibiting the precision determi- 
nation of the hadronic quantities, and of f B in particular. 

Lattice QCD is the closest to our ultimate goal, first- 
principle QCD calculation of the matrix element (|2jl. Al- 
though the calculations are based on the numerical simula- 
tions of the QCD vacuum fluctuations, the great bonus is 
that the only parameters appearing in the computations are 

2 See ref. [ 4 1 for the most recent update and for the full list of references. 
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those that appear in the QCD Lagrangian, namely the bare 
strong coupling and quark masses. A strategy to compute 
fs on a 4-dimensional lattice (L 3 x T) is very simple and 
can be summarized in 4 steps: 



350 



1. Generate an SU(3) gauge field configuration 1i (by , — i 300 
using the Monte Carlo technique); > 

^5 250 

2. For each time slice, t e [0, T), in the background <5 
field produced in 1, calculate the correlation func- ^ 
tion, ^^—-^ 
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2<0| eWroy 5 ?WAj(0)|o>. 
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Repeat steps 1 and 2, a number iVconf. of times, N C onf. 
refering to the number of independent gauge field 
configurations. From the average over the statistical 
sample, extract /# as 



<^A (x)Aj(0))« =° 



l<Q|<2y r 5 g|//>l 2 

2niH 
l -f>m H e- m »< + 



-m H t 



(4) 



Figure 1. Evolution of the lattice estimates of f Ds , as obtained by 
various lattice groups. List of papers corresponding to the above 
results, organized chronologically, can be found in refs. 1|7||8|. 
The last result, from ref. [|5|, is the final lattice value obtained in 
the quenched approximation. 



(iii) the simulations are made at several (four) small lattice 
spacings "a", which allows for a smooth extrapolation to 
the continuum limit (a — > 0). The final result of ref. [ 9 1 is 



where the (euclidean) time is large enough for the 
higher (heavier) excited states to be indeed sup- 
pressed w.r.t. the lowest lying one; 

4. Repeat 2 and 3 for several different light and heavy 
quark masses, m q and mg respectively. 

This simple procedure is very demanding in practice: to 
avoid large lattice artefacts, while working with reason- 
ably light quarks, the size of the lattice box (L) must be 
very large. Simultaneously, and to resolve the propaga- 
tion of the heavy quark, a very small lattice spacing ("a") 
is needed. Compared with the physical quark masses, the 
currently available computers allow us to work with 



f Ds =252 + 9 MeV. 



(6) 



< niQ < mi, , and ffl B /j < m q <m s 



(5) 



In other words, we can directly compute only the fo s decay 
constant. 

In fig. [2 we present the evolution of the results obtained 
from the lattice QCD simulations over the past 15 years. 
All results are obtained in the quenched approximation. 
The last value is actually the final quenched lattice esti- 
mate for this quantity. It is obtained after implementing 
all the important elements of theoretical progress made in 
the 90's, namely: (i) the quark action and the axial current 
are improved, providing a better scaling to the continuum 
limit; (ii) the local axial current defined on the lattice is 
matched to its continuum counterpart non-perturbatively; 



The fact that the Z?-quark cannot be resolved on the avail- 
able lattices opened three options: 

(a) Compute for the accessible mq and then from the 
fit of $}{mo) = C(m.Q)fu s/itih, in l/m#, extrapo- 
late to the B-meson mass. <!>(/«£>) is the quantity 
that scales with the inverse quark(meson) mass as 
a constant. The sizable 1/mg corrections are de- 
termined from the lattice data. The resulting val- 
ues, however, have large errors, mainly due to un- 
certainty of whether or not one includes the terms of 
0(1 /m"^ 2 ) in the extrapolation. Besides, the larger 
quark masses may induce large lattice cut-off arte- 
facts that are hard to quantify. 

(b) Difficulties in controlling the systematic errors of the 
strategy (a) incited many groups to attempt treating 
the heavy quark in an effective theory. Computa- 
tion of the matrix element in the static limit of 
HQET I 10 1 turned out to be very difficult, mainly 
because of the very poor signal [ II II . That prob- 
lem (of the poor signal) was circumvented by includ- 
ing the 0(1 //hq) corrections in both the Lagrangian 
and in the operator (axial current) by discretizing 
the NRQCD [[H). In spite of its benefits (working 
with large quark masses) this method cannot be used 
for the precision computation since, on the lattice, 
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the 1/niQ terms become <9[l/(amg)], where "a" is 
the small lattice spacing. In other words, the con- 
tinuum limit does not exist. Besides, in both ap- 
proaches (HQET and NRQCD) the non-perturbative 
renormalization of the axial current was not fea- 
sible. The third effective approach has been de- 
veloped by the Fermilab group \ 1131 . It consists 
of pushing the propagating heavy quark, Q, over 
the lattice cut-off and then expanding in powers of 
1/m.Q [not l/(amg)]. The key in that procedure is 
to match the relativistic with non-relativistic energy- 
momentum dispersion relations, where the mismatch 
is accounted in distinguishing the masses that ap- 
pear in the non-relativistic expansion as M- stat i C ", 
M-kinetic", and so on. That matching is typically 
made non-perturbatively while the renormalization 
is made only perturb atively. 3 

(c) Combine the value of <t>(yi) obtained in the static 
limit (1 /mg = 0), with the decay constants computed 
with the directly accessible heavy-light mesons, and 
interpolate to l/m/,. The main obstacles in follow- 
ing this strategy, as mentioned above, are the poor 
statistical quality of the static correlation functions 
and the missing non-perturbative evaluation of the 
renormalization constant. 

The overall agreement among results for fs,, as obtained 
by using various strategies, was quite impressive. How- 
ever, none of the approaches was fully satisfactory to pro- 
vide a precision measurement (even within the quenched 
approximation). This is why the errors remained essen- 
tially unchanged over the last 7 years, as can be seen in 
fig. 13 where we plot the evolution of the world average of 
the quenched lattice estimates of fg W151 . 



2 New ways to get f B (quenched) 
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Figure 2. Evolution of the world average of the (quenched) lattice 
estimates of f Bs , reported at the annual "Lattice" meetings r 1 1 51 . 
In the average enter the results obtained by all the strategies to 
treat the heavy quark on the lattice that are discussed in the text. 

2-2n( fl ) — * ^2«(0) = c 2 „. Throughout the calculation, they 
set T = 2L, and extract the decay constant precisely at T/2. 

The actual calculation with this clever procedure was re- 
cently completed in ref. [ 1171 . Let us first focus on fo s . In 
the first step, they work with Lq = 0.4 fm, so that at t — Lq 
they certainly cannot extract /d (Lo), but rather 



/o ff (^o) 



1 + 



2\ 

n 1 



m D,fo, 



-0n D (n)-mD,)Lo 



(7) 



where the sum runs over all excitations that are not heavy 
enough to be exponentially suppressed at (small) Lq. After 
doubling the size, Lq — > 2Lo, the exponential suppression 
will be more efficient, so that after 2 or 3 steps all excita- 
tions will die out. After performing the continuum extrap- 
olation in each step, they obtain 



The two results, reported this year, actually supersede all 
the previous quenched calculations. 

2.1 Method of finite volumes 

To simulate the b-quark directly on the lattice, one can ac- 
tually set a very small lattice spacing, but then the physical 
volume becomes too small to accommodate the propaga- 
tion of the light quark. To get around this problem, the 
"Tor Vergata" group [ 16 1 proceed as follows: (i) in a small 
box of side L = Lq, they compute /#(Lo); (ii) they double 
the side of the box, at fixed lattice spacing a, recompute 
f H (2L ), and take the ratio, Z 2 (a) = f H (2L ) / f H (L ); (iii) 
step (ii) is then repeated n times until 1o. n (a) — > 1; (iv) the 
whole procedure is redone for several values of the lattice 
spacing a, allowing the smooth continuum extrapolation, 



ff (L ) = 644 + 3 MeV, 



0"2 = 



cr 4 



/g ff (2£o) 

ff s (Lo) 
ffi(4Io) 



= 0.414(3), 



0.90(2), 



3 For details and a complete list of references, please see refs. 1 1141 . 



f D , = ff(Lo)o- 2 o-4 = 240(5)(5) MeV , (8) 



thus in good agreement with the value obtained in ref. 
The main reason why cr 2 is so different from 1 seems to be 
the presence of the non-decoupled excited states. This can 
be seen from what the authors of ref. [ 11711 present: they 
show that S2„(a) depends only weakly on both the heavy 
and light quark masses, and they verify only a tiny depen- 
dence on the lattice spacing. 

When working with the B s -meson, /J ff (Lq) can be com- 
puted directly, but doubling the side of the box to com- 
pute cr 2 „ becomes unfeasible. Instead, one can afford to 
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compute cr-inimQ), for m c < wiq < nib 12, and then reach 
0"2n('«2)X through an extrapolation. For that purpose they 
assume the heavy quark expansion in the finite box, so that 



(9) 



where the constant (3 naively scales as 1/(2" 1 Lq). Their 
data indeed verify eq. leading to 



/™-(Lo) =475 + 2 MeV, o- 2 (m b ) = 0.417(3), 

<r A (m b ) = 0.97(3), 
=> f B , = f^ (Lo)cr 2 (m h )cr 4 (m b ) = 192(6)(4) MeV . (10) 

The second error indicates the combined systematic uncer- 
tainty due to renormalization constants and to the contin- 
uum extrapolations. 

2.2 Combining with the static f B , 

As we mentioned before, the value of the decay constant 
in the static limit was plagued by two main difficulties: 1. 
the renormalization of the axial current in HQET on the 
lattice; 2. bad signal-to-noise ratio. Both problems have 
been solved recently. 

By a judicious choice of the renormalization condition, the 
authors of ref. [ fTSl provided a solution to the first prob- 
lem (see around eq.(2.15) of that paper). They devised the 
method of subtracting the power-divergent residual mass 
counterterm non-perturbatively, which then allowed them 
to apply the standard non-perturbative renormalization pro- 
cedure in the Schrodinger functional scheme [ 1191 . 

The solution to the second problem came recently too [ 20 1 . 
In the static limit, eq. (0 becomes 

c cont. (f) = ^i^^S^^ig^drA^A 
' l-l \ 

-^5,(d,()f]f/ (T)>, (id 



C latt (f) = <Tr 



r=I 



where S q (0, f) is the light quark propagator at x — 0, while 
the Wilson line on the lattice is simply the time-ordered 
product of link variables. Empirically, the statistical quali- 
ty of C latt (f) becomes much better if some kind of "fatten- 
ing" of the link variable is made: Uq(t) — > t/g at (T). 4 

After extrapolating to the continuum limit and converting 
the static HQET result to QCD, from ref. [ 1201 we learn 
that, 



/|' at = 225 + 16 MeV. (12) 

4 Fat link is obtained by averaging over the "staples" made of links that 
are the first neighbours to a link which is being fattened. Among various 
recipes to do that, the most efficient appears to be the one proposed in 
ref. l ETl . 



That result has then been combined with those that are ac- 
cessible directly, with the propagating heavy quark of mass 
(fling) < 0.65 (in lattice units) as to avoid larger discretiza- 
tion errors. After extrapolating to the continuum at each of 
the simulated heavy quarks, they show a very smooth linear 
behaviour of the combined set of data. Their preliminary 
result, presented in ref. f 1221 . is: 



f B , = 206 ± 10 MeV. 



(13) 



The bonus of that result is obviously a much better pre- 
cision, since the extrapolation to b-quark is now replaced 
by the interpolation. Notice however that the slope (i.e. 
1/mg-dependence) of <E>g = C/h, ym^, quoted in ref. [ 
1221 . is completely compatible with the previous calcula- 
tions in which the strategy with propagating heavy quark 
was used at a fixed value of the lattice spacing. 

2.3 Status and perspectives of quenched fo s and fs, 

After 15 years of computing the heavy-light decay con- 
stants on the lattice, in the quenched approximation, we 
are finally in a position to quote the precise values. From 
the results discussed in this section, I conclude 



fi 



■N t =0 



245 + 6 MeV . 



f N r =0 

Jb, 



194 + 6 MeV. (14) 



It should be stressed that these values are obtained by fix- 
ing the lattice spacing to ro = 0.5 fm, which, however, is 
an assumption. Using other quantities, such as fx, m.K>, 
leads to different values, which amount to a systematic er- 
ror of about 5%. The above results are obtained by work- 
ing with fine-grained lattices, by implementing the non- 
perturbative renormalization and after a smooth extrapola- 
tion to the continuum limit. Alpha group plans to attempt 
a computation of the l/m& corrections to the /^ tat , to check 
whether these are enough to feed the linear dependence that 
they observed after having combined the static and the re- 
sults obtained with relativistic (propagating) heavy quarks [ 

ED 

The most challenging issue that remains to be studied and 
understood is to assess the error induced by the quenched 
approximation. The knowledge gained about that system- 
atics is what we discuss next. 

3 Unquenching f Ds and f Bs 

To unquench the lattice results means that one has to in- 
clude the effects of dynamical quarks (those popping up 
through loops in the background gauge field). To incorpo- 
rate the light dynamical quarks is, in principle, possible, 
but very costly in practice. The most tested and widely 
used robust algorithm, Hybrid Monte Carlo (HMC) r !251 . 
allows one to include two light degenerate quarks of the 
mass close to the strange quark (please see [ 24 1). Getting 
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to the quarks lighter than half of the strange quark mass is 
nowadays impossible, if we use the standard Wilson quark 
action and keep the finite volume effects under control. To 
get over that limit, as well as to include the third (strange) 
dynamical quark, a substantial progress in algorithm build- 
ing is badly needed. An important step in that direction has 
been made in ref. I T25I . 

Up to now, the partially (un)quenched computations are 
made by using the Wilson quark action for the light and 
one of the effective approaches to treat the heavy quark. 

By confronting the results of quenched and unquenched 
studies (with Nf = 2 degenerate sea quarks) obtained by 
using the NRQCD treatment of the heavy quark on the lat- 
tices with fl"'»2 GeV [|26]|27]|28), the following effect 
of the dynamical quarks can be deduced : 



CP-PACS: 



JLQCD: 



f Nt= 
JB, 



1.10(5), 



f 

JB 



Z„* 1.13(5). 



JV f =0 



(15) 



With the Fermilab treatment of the heavy quark, instead, 
one has rl291P0l: 5 



f N f =2 



C N,=2 



CP " PACS: %=o =1-14(5), ^ = 1.07(5), 

Jb, jd s 

f N f= 2 f wf=2 

MILC: Jl—~ i.Q9(S) t ~ 1.08(5) . (16) 



f- 



JV f =0 



f Nf=0 



The average of the above results, combined with the values 
given in eq. 14l >. leads to 



f N f =2 



'— = 1.12(5) => f"<- 2 = 217 ± 12 MeV, 

/-JVf=0 J fij 
B, 

JD. 



r 

JB, 
pN f=2 



- = 1.08(5) => f^- 2 = 265 ± 14 MeV. (17) 



JV f = 



JD 



Unaccounted for is ~ 5% of systematic uncertainty, which 
we should keep in mind, due to the scale setting when con- 
verting from the lattice to the physical units. A detailed 
study of that uncertainty in the unquenched (Nf = 2) simu- 
lations has not been made so far. 

Finally, we are interested in the situation in which Nf = 
2 + 1, where the extra flavour would correspond to the 
strange sea quark. With Wilson fermions such a study is 
extremely expensive, although the algorithms to do such 
simulations already exist f 13 II . Instead, a study with the 

5 The MILC group made the most extensive unquenched (Nf = 2) study 
of the heavy-light decay constants. Please see ref. 1 for the detailed 
report on their findings. 



so-called staggered light quarks has been made [ 32 1. Al- 
though the method is relatively cheap to implement, one 
of the unsatisfactory features can be formulated as follows. 
With the staggered action, each dynamical quark flavour on 
the lattice comes in four "tastes" (copies), and it is unclear 
how one can relate such Dirac determinant A sta ggered, to the 

/ \l/4 

desired one. A current practice of taking + 1 A sta gg er ed) 
induces non-localities that are potentially problematic, and 
the whole formulation may not correspond to QCD. The 
proponents of the method will hopefully study this issue 
more carefully. Having that off my chest, I can now com- 
pare the quenched and unquenched results of ref. [ 1321 . in 
which the NRQCD treatment of the heavy quark has been 
used, to quote 



(/*= 2+1 //£ f=0 )*1.15, 
i.e. fully consistent with the result in eq. H7\ . 



(18) 



In conclusion, the current values of the pseudoscalar decay 
constants, with the light quark s, are: 



fa. 


= 217 ±12 ±11 MeV, 




fo, 


= 265 ±14 ±13 MeV. 


(19) 



4 fa J 'h m and f Ds /f Du/d 

So far we discussed the case with the strange light quark 
in the heavy-light meson. In particle physics phenomeno- 
logy, however, more needed is the information about fa uld , 
fD uld , and about the SU(3) breaking ratios, /b,//b„ w and 
fo. If D m ■ The best known example of how important is the 
knowledge of these quantities is the one given in eq. Q. 

As we already pointed out in the previous sections, the 
physical m M /j is too small to be simulated directly on 
the lattice. Currently feasible are the light quarks whose 
mass, with respect to the physical strange quark mass 



(r = m 9 /mf ys '), lie between 



1/2 < r <2, 



(20) 



implying the necessity for quite a long extrapolation to the 
physically relevant limit, r — > r u /4 - 0.04 \ 1331 . This 
is where the staggered quark action has a great advantage 
over the standard (Wilson) one: due to the explicit chiral- 
ity, with staggered quarks one can reach quarks as light as 
r ~ 1/4 [ 32 1, implying a better control over the extrapola- 
tion to r u/d . 

Within the accessible range of the light quark masses i2Q\ . 
the results of both quenched and unquenched simulations 
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(with Nf = 2) suggest a pronounced linear dependence on 
the light quark mass, 



^ = 1 + X h -(l-r), 

JB q 



(21) 



with little or no room for a term oc r 2 . For the quenched 
value of the slope Xj,, I will take the one quoted in ref. r i!7l . 
namely X^'~° = 0.13(2)(1). To get the unquenched slope, 
we first compile the available results for the SU(3) break- 
ing ratio as computed with iVf = 2, divided by its quenched 
value, where both results are being obtained from an ex- 
trapolation of the form (12 1 1 . Those results are plotted in 
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A common practice "to use a ruler" and extrapolate the 
observed behaviour to the physical r — > r u /d is, however, 
dangerous because one is getting deeply into the region 
dominated by the spontaneous chiral symmetry breaking 
effects, so that the observed linear behaviour of the decay 
constants w.r.t. the change of the light quark mass may 
be significantly modified. That problem was first pointed 
out in ref. \ 1351 . and recently in \ 1361 . It is therefore 
of paramount importance to do the computation with the 
light quark as light as possible, and check for the deviation 
from the form (12 1 i . In doing so, the finite volume effects 
should, of course, be kept under control. That is where 
the e-expansion should be applied in the way similar to 
what has been done recently in ref. f 1371 . 6 Since such an 
unquenched study is not around the corner, we may rely 
on the chiral perturbation theory (ChPT) to guide the (chi- 
ral) extrapolations. By using the Lagrangian in which the 
ChPT is combined with the static limit of HQET, the chiral 
logarithmic corrections to the SU(3) breaking ratio of the 
decay constants has been computed in [ 38 1 



R 



ChPT 



1 + Ji^( 3 "'^A, 



-2,4 log/ 



m 2 ^ 

2 , 1 

m n \og — 



+ 75— ( m K 



ml), 



(24) 



where K stands for the unknown low energy constant. 



Rg d is independent of the chiral symmetry breaking scale 
A x 1 GeV. This formula can be cast into the form ready 
for extrapolation by using the GMOR and Gell-Mann- 
Okubo formulae, 



m n = 2B m s r , 



2B m s 



r+l 



m — 2B(,m s - 



r + 2 



Figure 3. Ratio of the slopes which appear in eq. 12 1 i in light 
quark mass, as obtained from the quenched and unquenched sim- 
ulations, by using the NRQCD and the Fermilab treatment of 
the fo-quark. From up to down, the results correspond to refs. [ 
EHSQ.[|2D> [H3 and [|29). 



fig. The average is 



N{=2 



(fBjf Bd y 

(fBjfB d ) N '=° 



= 1.02 ±0.03, 



or, the SU(3) breaking ratio is only slightly sensitive to the 
switch from Nf - — > 2, for the light quark mass in the 
range indicated in eq. J20I . Notice also that a preliminary 
study with staggered fermions indicates that this feature 
persists even after going to Nf = 2 + 1 . From the above 
discussion we can now easily extract the slope: 



X?'- 2 = 0.15(4). 



(23) 



After noticing that 2B m s = 2m 2 - m 2 K = (684 MeV) 2 , 
eq. ( I24> can be written in terms of r, which we will then 
call /?£ hPT (r). Notice that the last term is oc (1 - r), much 
like in eq. (12 It in which the slope was obtained from the 
fit to the lattice data, and which we will refer to [eq. ( 12 It . 
that is] as Rg(r). The main trouble actually lies in the fact 
that the chiral logarithms are multiplied by the large coef- 
ficient (1 + 3g 2 ) ~ 1.8, 7 so that the application of the form 



R)f {r) to extrapolate to r u /a results in a large shift, com- 
(22) pared to what we obtain from the naive extrapolation, i.e. 

from the result of the application of /?^ tt (r„/d). Since we do 
not directly observe the chiral logarithms on the lattice, the 



6 It would be very interesting to actually verify the predictions of quenched 
ChPT, which trouble the lattice QCD community so much. Using the 
method of ref. [ 1371 . one could check whether or not the (divergent) 
quenched chiral log term in the fitlfn coincides with the QChPT pre- 



1 + r. 



■ log r, where mo a 0.65 GeV, is the mass of the 



diction, — 

6(4tt/) 2 1 - 

light quenched T]' state. 

1 g is a coupling of the Goldstone boson to a doublet of the heavy-light 
mesons (e.g. gg»s T ), whose determination will be discussed later on. 
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formula R^ hVT (r) can be used only below the region of the 
quark masses covered by our lattice simulations, namely 
for r < tm - 0.5. It can also be argued that R^ hPT {r) is not 
appropriate for the 'pions' as heavy as m-„" = 0.48 GeV 
(o r M = 0.5), and the form ^ att (r) 

can be extrapolated to 
lower rM, before including the ChPT formula in the extrap- 
olation. Clearly, the point ru, at which the smooth match- 
ing can be made, 



R Bslq (r) 



&(r-r M )R l *l(r) 



+§{r M - r) 

3 ChPT 



X b j (r - r M ) 



(25) 



will result in different shifts. By taking = 0.75, a pure 
SU(3) breaking effect, / Bj - 1 , gets enhanced by 100% [ 
l36l. If, instead, we take r M = 0.25, the shift is still +30%. 
This is illustrated in fig. |4] An equivalent way to see that 
feature has been discussed in ref. f 1391 . where the chiral 
loop integrals are computed by introducing the hard cut- 
off. The uncertainty on the value of the cut-off scale cor- 
responds precisely to the uncertainty of the position in tm 
discussed above. To get around that uncertainty, before we 



Linear extrapolation 

Linear+ChPT (r M =0.25) 

Linear+ChPT (r =0.75) 




Figure 4. Illustration of the shift induced by the chiral logarithms 
in the extrapolation from the light quark masses accessed from the 
lattice, r e (0.5, 1], to the physical point, r u j c \ = 0.04. The sensi- 
tivity of the resulting f Bs jf Bd on the choice of the matching point 
r M [see eq. <25H is also shown. The phenomenon is illustrated by 
using X b = 0.12(2)(3), and g = 0.5. 



are able to work with the light quarks sufficiently close to 
the chiral limit on the lattice, in ref. [|40J it was proposed 



to study the double ratio 



(26) 



There are two reasons why this is advantageous: 

A. The coefficients multiplying the chiral logarithms in 



flChPT s f* = j + 



fn 



4(4*/) 

2 



5m l\ g-* -2m \\ og 



A, 



A. 



2 m n 
-3m, log — 



8L 5 (A X ) 

f 



{m\-ml), (27) 



have the same sign and are almost of the same size 
as those appearing in eq. ( 1241 . 8 Therefore, in the 
ratio, the logarithms will cancel to a large extent and 
the form 



R(r) = A + B ■(! - r) 



(28) 



can be used in extrapolations, with a very small error 
due to the uncancelled logarithms. Using an equa- 
tion similar to ( I25> . but for R(r u /d), we notice a weak 
(negligible) sensitivity on the choice of r^, because 
the large logarithms are cancelled. 

B. The physical result is obtained after multiplying 
R( r u/d) by fx/fji, which is known from experiment, 
namely 1.22(1) [HQ. 

Concerning /f„, it is important to stress that, compared 
with the experimental value, the quenched lattice estimates 
obtained after extrapolating linearly, 



R 



latt 



= 1 + X„-(l-r), 



(29) 



always lead to a small value. For example, from the 
SPQcdR data [ 42 1, in the continuum limit, I get Xg' =0 = 
0.12(2), which then amounts to, fx/fn - 1 = 0.11(2), i.e. 
50% smaller than the experimental value. The folkloric ex- 
planation for this "spectacular failure" was/is the use of the 
quenched approximation. JLQCD collaboration recently 
showed that this is, however, not totally true. From their 
precision quenched and unquenched (iVf = 2) computa- 
tion, the results of which are reported in ref. r 1431 . I read 
off 9 



X 



Nf=0 



0.125(7), 



0.154 (30) 



This statement obviously depends on the value of the coupling g, which - 
in the fo-quark sector- 1 take to be gt = ().52(7)(7), which will be discussed 
in the last part of this review. 

'JLQCD also made a preliminary study with Nf = 3, indicating that 



M=3 



»JVf=2 



8 
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The missing piece that would help getting to (//f// ff ) exp ' 
is likely to be the one due to the chiral logarithms. This 
situation is precisely the opposite to the one that incited 
the ChPT practitioners to make the consistent NLO calcu- 
lations, because the chiral logs alone give fxlfn — 1.14, 
whereas at NLO the low energy constant L5 [see ea. (!27H . 
is fixed so as to reproduce (f K l '/ 7r ) exp ' [El . 10 On the lat- 
tice, instead, we see the explicit linear quark mass depen- 
dence but there is still no clear evidence for the presence of 
the chiral logs (most probably because we are not working 
with sufficiently light quarks). To get the lattice estimate 
of /kZ/tt, we then face the same problem as before: chiral 
logs are large and the results of extrapolation by using an 
expression similar to eq. i25i would be strongly dependent 
on the choice of the matching point r^- 



5 5° - B° q mixing and £ 



The double ratio i26\ avoids all those headaches, and can 
be computed on the lattice directly. By combining X q = 
0.154 (+}*), with X b = 0.15(4), I get 



hi 



f exp 
■>K 

f exp 
.In 

1.21 



1 + X b ■ r u/d 



1 + X„ 



r u /d 



x "chi - logs'' 



: 0.05 ± 0.01 



(31) 



where the last error is due to the variation of the coupling 
g = 0.52(7)(7). A similar proposal, to consider yfrnsfB/fn, 
has been made in ref. f 1301 . However, one has to assume 
that the terms of 0(1 1 nib), which are known to be large, do 
not induce any extra chiral log dependence. Large 0(1 /nib) 
correction cancel in the /bJ/b ratio. 

A completely analogous procedure to the one sketched 
above applies to the charm sector. I will use X c - 
0.11(4), which is obtained from X* =0 = 0.09(1)(1) [ 
17 1, and the observation made in refs. r 1291 1301 . namely 
(fDjfD) N ' =2 /(fD,/fD) N ' =0 = 1-02(3). Together with the ex- 
perimentally measured g c = 0.61(6), and from the formula 
analogous to eq. d3 1 1 . 1 arrive at /dJ/d = 1.22(4). 

In summary, 



- 1.22 + 0.04, 

JD 

ifi = 1.21 ±0.05 ±0.01 , 
Jb 



(32) 



or, within the error bars, fajfa ~ /d J/d ~ hlfn- 



l0 As a side remark, one should stress that neither in ChPT nor in QCD 
sum rales the fic/fn is obtained directly, but rather after fixing the extra- 
parameters, i.e. not only changing the quark mass m„ j <-» m s (for the 
recent QCDSR analysis of ffclfn, please see ref. 1 1451 ). This is why 
a thorough lattice study, allowing a deeper understanding of the SU(3) 
breaking mechanism, is badly needed. 



The bag parameter that parameterizes the Z?" - B q mixing 



amplitude is defined by 



(B° q \b%(l -yd? ^(1-75 )q j \B q ) = ^I/IS , (33) 

i.e. normalized to the vacuum saturation approximation 
(Bvsa — 1)- This amplitude is very difficult to calculate 
by QCD sum rules because: (i) the calculation of the NLO 
corrections to the perturbative part of the spectral function, 
as well as to the Wilson coefficients multiplying the con- 
densate contributions, is quite involved; (ii) the stability of 
the sum rule under the variation of the threshold and the 
Borel parameters is very hard to achieve. Some progress 
concerning the first part of the problem was made recently 
in ref. f 1461 . So far, the lattice QCD studies of the ma- 
trix element fl33i are made by using the Wilson quark ac- 
tion. The unpleasant feature of the Wilson quarks is the 
lack of explicit chirality, which induces complications in 
renormalization, namely it allows the operator in eq. J33i 
to mix with other parity-even, AB = 2, dimension-six op- 
erators. The extra mixing is the lattice artefact and should 
be subtracted, preferably non-perturbatively. To a discus- 
sion provided in the Yellow Book from the previous CKM- 
workshop [ 2 1, 1 would like to add the recent proposal made 
in ref. [ 1471 . to combine the overlap light quark action 
(which preserves the chirality) with HQET, which enor- 
mously simplifies the renormalization procedure. 

In fig.|5l we show the results for the SU(3) breaking ratio of 
the bag parameter, as well as the absolute value of BsOiib) 
in the MS(NDR) scheme. Let us briefly discuss some im- 
portant points. In ref. [ 50 1, the static HQET result has 
been combined with the ones obtained with the propagat- 
ing heavy quarks. That was an important progress because, 
instead of extrapolating, one could interpolate to get the de- 
sired result. In doing so the matching of the QCD results 
with the static (HQET) ones has been made at NLO in per- 
turbation theory, which is important for two reasons:(l) the 
anomalous dimensions and mixing patterns are different in 
the full and effective theories (QCD and HQET), and to 
use the heavy quark scaling laws the matching to HQET is 
required; (2) to cancel the scale and scheme dependence 
against the perturbatively computed Wilson coefficient [ 
1511 [i.e. Cfi(ju) in eq. JD], one has to provide Ba sd (p), com- 
puted in the same MS(NDR) renormalization scheme. All 
schemes are equal at the leading order, and to specify the 
renormalization scheme, one must go beyond, i.e. to NLO. 
All that is done consistently in ref. f 1501 . However, those 
results are quenched, they are obtained at a single value of 
the lattice spacing, and the matrix elements in HQET were 
renormalized only perturbatively. Important step in tam- 
ing the quenching effects has been made in refs. [ 1281 1521 . 
where the heavy quark is treated non-relativistically. From 
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propagating heavy 

A prop. heavy + static limit 

■ NRQCD 



0.97(12) 



B B = 1.37(14), B B JB B = 1.00(2), £ = 1.21(6) . (35) 



0.98(5) 



0.98Q 



1— 



0.93Q 



0.910 



' • 



0.99(2) 



1.02(3) 



1.02(5 



0.87(6) 



0.84(6) 



Bernard et al. 



APE 

2000 



UKQCD 

2001 



SPQR 

2002 



JLQCD 

2002 



Since there is still much room for the improvement of the 
bag parameter results, I should also mention the recent 
QCDSR result, B^{m h ) ~ 0.95 [|46|, as well as the one of 
the chiral quark model, B B JB Bd = 0.93(5) [ 5!. 11 

6 g coupling 



0.84G 1 ) 

I ■ 1 



The coupling of the lowest lying doublet of heavy-light 
JLQCD [unq. N,J?f sons O'f = I/ 2 ") t0 a charged soft pion, g H - H n, is de- 



0.9 



B 



1.0 



i.i 



0.8 0.9 1.0 
B B ( m J 



Figure 5. SU(3) breaking ratio of the bag parameters and 
B^(m b ), as obtained from the lattice studies by extrapolating 
from the accessible propagating heavy quarks, by constraining 
the extrapolation with the static HQET result, and by using the 
NRQCD treatment of the heavy quark. The plotted results are 
taken from refs. r|48||49||8l|5o||28l|52|. 



the comparison of the results obtained in simulations with 
Nf - and A^f = 2, they see no difference between Bg f ~° 

and/?' 2 

We should reiterate that the light quark in BbXh) is ac ~ 
cessed directly, whereas Bb^h) is reached through an ex- 
trapolation. ChPT suggests no deviation of the SU(3) 
breaking ratio from the linear form. From ref. [ 1381 we 
know that the chiral log term reads 



B Bd 



= 1 + 




2(4tt/> 



[m 2 log m 2 -m;; log m 2 ] , (34) 



which, contrary to the case of the decay constants, is very 
small; no significant shift in the chiral limit is therefore 
to be expected. That gives us confidence that the chiral 
extrapolations do not induce important systematic uncer- 
tainties, although that issue will be resolved iff we are able 
to actually do the lattice computation close to the chiral 
limit. By taking the simple average of the results given in 
fig.0 we obtain B*f{m b ) = 0.89(3), to which we add 10% 
of uncertainty due to possible remaining: quenching, chi- 
ral extrapolation, continuum extrapolation effects. To con- 
vert from MS to RGI form, with A^ D = 213(32) MeV [ 

53 1, the conversion factor in B Btl/s = c(mb)B Bdl ,( m h), is 
c m (4.8GeV) = 1.54(1). The average SU(3) breaking ratio 
of the bag parameters gives 1.00(1), the error of which we 
double for the same reasons mentioned above. Therefore 
we have 



fined as 



(H(p')Mq)\H*(p)) = g H . Hn (q ■ e(p)) . 



(36) 



This coupling appears to be the essential parameter in 
the chiral Lagrangian of the heavy-light meson systems, 
Xim = gTr[HHy IJ y 5 A' 2 ], with g oc g H , Hn . The constant 
gH-Hn is important in describing the D — > niv and B — > nlv 
decay form factors, since the residuum of the dominant 
form factor at the nearest pole (i.e. m 2 D , and m 2 B „ respec- 
tively) is given by fH>mH'gH'Hn/2- The value of this cou- 
pling has recently been measured by CLEO in the charm 
sector, go-Dn - 17. 9+0. 3± 1.9 [ 55 1, which we then convert 
into the value that is expected to scale with the heavy quark 
mass as a constant (up to 1/rriQ corrections and higher), 
namely 



gD'Dn - 



2 y/m D m D . 



fn 



g c , =0.59(7), 



(37) 



where the index "c" indicates that the heavy quark is 
charm; g D , on the other hand, cannot be measured directly 
because there is no available phase space for the pion emis- 
sion (or absorption). A summary of the predictions of this 
quantity as obtained by using various approaches can be 
found in ref. [ 56 1. Here I will focus on the recent develop- 
ment. 

The experimental finding by CLEO was surprising in that 
the value for gwDn was much larger than the one predicted 
by the light cone QCD sum rules (LCSR) f BTl . To cure 
that discrepancy, ref. [ 58 1 proposed to include the first ra- 
dial excitations in the hadronic sum, i.e. below the con- 
tinuum threshold in the double dispersion relation. As a 
result, the LCSR value for gu-Dit becomes much larger and 
more stable against the variation of the sum rule parame- 
ters. In addition, the sum rule for the D^'-meson decay 
constant remains unchanged, whereas the D — > n£v form 
factor gets corrected in such a way that the result of pole 
dominance becomes closer to the experimental value for 
this form factor. 



"From the plethora of predictions given in ref. 1541 . 1 chose to quote the 
one for Bg I /BB d , which is very stable under large variations of the model 
parameters. 
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0.2 0.4 0.6 0.! 



CLEO (2001/02) 

0.61(1)(6) 

Orsay (2002) 

0.67(8)(5) 

KRWY (1999) 

0.36(10) 

Orsay (2003) 

0.58(6)(10) 
KRWY (1999) 

0.27(9) 

Orsay (2003) 

0.48(3)01) 



tice in the static limit of HQET. 12 Orsay group obtains 
goo = 0.48(3)(10), where the bulk of systematic uncer- 
tainty stems from various ways of smearing the source op- 
erators (the ones that produce H*-*') f 1611 . Thus the result 
for the slope is very close to the one predicted by LCSR, 
y lcsr , although the absolute value obtained on the lattice and 
from LCSR do disagree. With the static result in hand, 
one can interpolate to the Z?-meson mass, which results 
in gi, = 0.58(6)(10). 13 The present situation of the esti- 
mates based on LCSR and quenched lattice simulations is 
shown in fig. HO Besides the obvious necessity to go beyond 
quenching and to better control the chiral extrapolations, it 
would be nice if other lattice groups produced results for 
this coupling, even in the quenched approximation. 
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Figure 6. Present situation of the estimates of the g-coupling 
between the lowest doublet of heavy-light mesons and the pion. 
Plotted are the results obtained by using the light cone QCD sum 
rules (LCSR) [|57), (quenched) lattice QCD [EHUD and the 
experimental value [[33). 
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The first lattice computation of this coupling was made last 
year [ 59 1. The findings of that reference can be summa- 
rized as follows: 

1. From the simulations in the quenched approxima- 
tion, it appears that the finite lattice spacing and fi- 
nite volume effects are small; 

2. The coupling is computed from the transition form 
factor between the light quarks in H and H*, via 
the axial current, with the heavy quark being only 
a spectator. The coupling g q is almost insensitive 
to the value of the heavy quark mass for the heavy 
masses around the charm quark (see fig. 5 of ref. [ 
ED); 

3. The resulting value is g c = 0.67(8)Cg), i.e. it is large 
and in good agreement with gT P . 

Concerning point 2, the observation that the slope in 

gQ =goo (1 + r/m Q ) (38) 

is y latt as 0, disagrees with the LCSR prediction, y lcsr as 
0.6 GeV. It is important to say that LCSR came to that con- 
clusion by computing the couplings g b and g c , while on the 
lattice only the masses close to the charm quark are stud- 
ied. Since gb is highly important [see e.g. eq. J24H . it is 
necessary to have a better control of the slope y. To get 
around that problem, the Orsay group very recently fol- 
lowed the strategy of ref. [ 60 1 and computed g on the lat- 
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